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Abstract 

The time dependent Schrodinger equation with a time dependent delta function 
potential is solved exactly for many special cases. In all other cases the problem can be 
reduced to an integral equation of the Volterra type. It is shown that by knowing the wave 
function at the origin, one may derive the wave function everywhere. Thus, the problem 
is reduced from a PDE in two variables to an integral equation in one. These results are 
used to compare adiabatic versus sudden changes in the potential. It is shown that 
adiabatic changes in the potential lead to conservation of the normalization of the 
probability density. 
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Introduction 


Very few cases of the time dependent Schrodinger wave equation with a time 
dependent potential can be solved exactly. The cases that are known include the time 
dependent harmonic oscillator ,[1-3] an example of an infinite potential well with a 
moving boundary ,[4-6] and various other special cases. There remain, however, a large 
number of time dependent problems that, at least in principal, can be solved exactly. 

One case that has been looked at over the years is the delta function potential . The 
author first investigated this problem in the mid 1980s [7] as an unpublished work. This 
was later included as part of a PhD thesis dissertation [8] but not published in a journal. 
In that work the Schrodinger wave equation with a strength that varies as a function of 
time was investigated. That work was more extensive than what we present here in that it 
also included the periodic time dependence and a delta function array with time 
dependence. In the time that has transpired since then, much (published) work has been 
done in relation to time evolution delta function potential problems. Most of that has 
focused on the time evolution/propagator derivation of the case of one or more delta 
functions with constant amplitude in time. In contrast, this work deals mainly with a 
potential with a strength that varies in time. 

There are two types of time dependent problems that are commonly investigated. The 
problem that has received the most attention over the years is the simple scattering 
problem where one seeks what amounts to a steady state solution where the potential has 
a periodic dependence on time. In the case of the delta function potential, this problem 
was investigated where the strength varied sinusoidally in time[9]. Using Floquet 
formalism, they were able to derive a transmission coefficient as a function of driving 
frequency and amplitude. Other examples of scattering problems involving delta function 
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potentials also exist in the literature [10]. The actual solutions, however, are simple plane 
waves. 

The other time dependent problem is where one starts with an initial state and allows 
that to evolve in time. These are the diffusive solutions to the Schrodinger wave equation. 
This and the closely related problem of finding the propagator have been investigated in 
some depth in relation to the one or more delta function potential [11-16]. Out of those 
cases cited only one of them involves a potential with true time dependence [13] and that 
is the case of two delta functions that move apart from one another at constant velocity. 
There are two common approaches one may use. One method is using a path integral of 
the Feynman type and the other is using Laplace transforms. We use the latter though the 
two methods are closely related. 

Formulation of the Problem 

We would like to solve the Schrodinger equation with a time dependent potential , 

"V^ + 2 c(t)8(x)y/=i\i/ t 
— > 0,jc — >+°o 

By integrating Equation 1 over the discontinuity at x=0, we find that the problem 
reduces to, 


r 

yr (0V)-yr J[ ((r,f) = 2c(f)yf(0,f) (2) 

yr -> 0,;c -> ±°o 

This is just the case of zero potential with an additional time dependent boundary 
condition. If we restrict ourselves to t>0, we may take the Laplace transform of Equation 
2 


-V xx = is\l/-iV{x, 0) 

yr x (0 + ,s) - y/ x (0~,s) = 2L[c(f)y<0,f)] 
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We now solve Equation 3. However, we must consider the homogeneous solution in 


order to satisfy the cusp condition at the origin. The solution is 

i f/(x,s) = a(5 , )exp^'V^|xjj + | dx' exp(i-\[is\x - 

where a(s) is an arbitrary function of s, yet to be determined. By applying the cusp 
condition at the origin to Equation 2.4 we find 

\ff x (0 + ,s) - y/ x (0~,s\ = 2 i-\[isa(s) = 2L[c[t)y/[0,t)] 

L[c(t)\if(0,t)] 


a(s) = 


i'Jis 


(4) 


(5) 


We now set x=0 in Equation 4 and use the results of Equation 5 to find 

y/(0,s) = t ^ ^ + — 1 J dx'exp(i'^is\x'f\\i/[x' ,0) (6) 

His 2V is ' ’ 

Using the convolution theorem for Laplace transforms [17], we invert Equation 6 to 


find 


<//((),;)- = - 1 J dx'exp 

i4in o Vf-f 2 -Jim 


J dx'exp 

r 


4 1 

K 7 


Mx',0) 


(7) 


Equation 7 is a Volterra integral equation that determines the wave function at the 
origin. Once the wave function at the origin is known, one may find the wave function 
everywhere by inverting Equation 4 which has the result 


¥{xd) — . r — j 
l\ IK Q 








2-yJint 


dx'ex p 


.{x-xf 
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Mx', 0 ) ( 8 ) 


Equations 7 and 8 completely determine the problem. As we can see knowing the wave 
function at the origin is equivalent to knowing the wave function everywhere. 
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Some Specific Cases 

Derivation of the Propagator for Constant c and Jumps 

In some situations it is advantageous to represent the solution in the form, 


= J dx'G(x,x' ,tyf/[x' ,0) (9) 

where G is the propagator. To find the propagator, we use Equations 5 and 6 to find, 


| d[x'exp(z>/fs|x'|)i/r(x',0) 

l\ IS ^^Sy/Si^lC 


( 10 ) 


By Equations 4 and 10, we find, 


C ex p( 


X + x' 

)) expl 

(i~Jis\x - x'j 

2 + ‘\ficj 


2-\[is 


( 11 ) 

By making the substitution k =-\lis we see this is the same propagator found using the 
Feynman approach [11], We may now invert Equation 1 1 to obtain, 


G(x,x',t) = -^exp[c(|x| + |x'|) + zc 2 t]erfc 


1 n~ \x\ + \x' ^ 

cslit + 


V 


2-\[it 


1 


l4int 


exp 


( ( ^ 

(x-x ) 

i- — 

At 


( 12 ) 


Non-adiabatic Jumps 

Although this example is for c=constant, we may use this to generate a solution for a 

sudden, non-adiabatic jump in the potential. This is quite simple to do. If we choose c=- 
c 0 , where c 0 is a positive real number and take y/(x,0) = yfc^ exp(-c 0 |x|) as the initial 


state, the system evolves with a probability density that doesn't change in time except for 
phase. That is because we have chosen the bound state for the attractive potential case. 
One way to model a sudden jump is to chose an initial state of 
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y/(x,0) = exp(-c 0 |x|) , but let it evolve thereafter with c=-c', where c’>0. It can be 


shown that in this case [8], 

2 c ^ c o 


yr(x,t)^ — — exp(-c' x )exp(ic t),t— »°° 


c o + c 


The total probability of it radiating away is, 


(13) 


1-(| Vr(x,t)| 2 ) = 


( c 0 -c' 


y c o + c j 


,t ° 


(14) 


We see that sudden changes in potential lead to irreversible changes in the normalization 
of the probability density as we would expect. 


Adiabatic Limit 

Interestingly enough, we may also generate the adiabatic case as well using similar 
reasoning. The probability of it remaining in the bound state after one jump is 


(|yr(x,t)| 2 ) = 



,t — > °o 


(15) 


Lets suppose now we break it up into two steps. We first take a small half step, let it 
evolve for a long time, then take a second step and let it evolve again. The result is 


(kM 2 ) 


4CQC, 4c t c 2 
(cq+cJ 2 (c^cj 2 ’ 


t —> °° 


Here, c n = c 0 + n(c'~ c 0 ) / 2 . After N such steps we have, 


(16) 


<k(x,, w p)=4"n ; 

where c n = c 0 + n(c'~ c 0 ) / N . The result is 


4 c n c„ + i 

n=0 ( C n +C n+l) 




(17) 


6 



(|i/A(x,t,N)| 2 } = 

4 c 0 (c'+ c 0 (-1 + N))N r(a, - 1 + N) T(a 2 - 1 + N) ( r(a 3 ) 


(c'+c 0 (-l + 2A r )) 2 r(a t ) r(a 2 ) ^r(a 3 -l + iV) 


,t — > oo 


where 


a x = 2 + — j- 5 — A^, 


a 2 = 1 + — °—N, 


a =l + _^iv. 
2 c'-c ft 


By taking the limit as N approaches infinity we find 


lim(|vr(x,tJSf)| 2 ) = l, 


t — > oo . 


This implies the probability of the particle radiating away is 0. Thus, for adiabatic 
changes in the potential, the normalization of the probability density is conserved, as we 
would infer from the Bom-Fock theorem [18] on the conservation of quantum numbers. 


Solution for c Linear in Time 

We now consider a potential that has the following form, 

c{t) = -c 0 ,t< 0 
c(t) = -c 0 - at,t > 0 

Here, a is a real number. We take the bound state for c=-co as the initial condition, 
which in this case is y/(x,0) = yfc^ exp(-c 0 \x |) , where a and c 0 are real constants and 

c 0 >0. By Equation 6 and using L[ t f(t)]=-d (L [f(t))]/ds, we find, 


y/ s (0,s) - — (n/is + c 0 )yf(0,s) = 
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We solve Equation 22 with the condition that the Laplace transform vanishes at s=°° 
and that it must also be bounded at ±i°°. The last condition is necessary for the inversion 
integral to exist. The solution is, 


fr\ \ W C 0 ( 2 / . \3/2 Cq 

w(0,s) = — - — exp — (is) +—s 
YK ' a {3a y ’ a , 


] ds '~ 


2 / 


\ 2 / . ,\3/2 Cq 

exp (w ) — -s 

, 3« v y « J 

c 0 - i\[is' 

2 / . ,\3/2 C n 


,a >0 


(23) 


yr(0,5) = — 0 exp 
a 


3a 




^ 1 . |\3/2 <-Q 

exp ) — -s 

1 3 a y ' a J 

c 0 - i^fis' 


,a< 0 


This expression cannot be inverted in terms of elementary nor special functions. 
However, It can be inverted using the integral, 


1 - 

y/(0 ,t) = — J dsy/(0,s)exp(st) (24) 

2jti -i~ + p 

where P is a positive constant chosen to be greater than the right most pole in the 
complex plane. In this case it would be chosen to be greater than 0. 

One application for this would be as in the previous example where one could 
possibly compare adiabatic versus instantaneous change in the potential. This linear time 
dependence was studied in some detail by Dodonov[19] for the Coulomb potential. 
However, we have thus far been unable invert the Laplace transform represented by 
Equation 23 analytically, though numerically it is not a problem. 


Solution for c Inversely Proportional to Time 

Potentials that collapse in time could possibly be used to model nuclear decay. 
Another example might be a collapsing symmetric electric multipole in 1-D. 

We choose c of the form, c(t)=a/t, where a is a real constant. In this situation, we 
must be very careful. The Laplace transform of c(t) doesn't exist. However, the Laplace 


8 



transform of c(t) \|/(0,t) does exist if we choose an initial condition that vanishes at the 
origin. With this restriction in mind, we proceed with the problem. We first use the 
identity, 

L[-f(t),s]=]ds'L[f(t),s'] (25) 

* i 

By Equations 25 and 6 we find 


Let, 



m(^)= J<is'yr(0,£') 

s 


Then by Equations 26 and 27 we find 


(26) 


(27) 


u s + u = J dx'exp(i^\x'tyii/(x' ,0) 


i-yjis 2 -\lis 

After solving Equation 28 we find 


(28) 


exp| 


*') 


x' 

+ 2 ia 


By Equation 29 we find, 

vHp, 


2 -yjis . 


x' 

exp x'j 


x' 

+ 2 ia 


We may now invert Equation 30 to find, 


(29) 


(30) 





x' 


x' 

+ 2 ia 


,.(*f 

4 1 




(31) 


We now seek \|/(x,t). To find this quantity, we first write, 


°a 

1 f/[x,t) = _ dx'G(x,x' ,t)if/(x' ,0) 


(32) 
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To find G, we use the results of Equations 30, 27, 5, and 4 to find, 


1 


G(x,x' ,s) = — j= 
2^ is 


2 ia 


|jq + 2 ia 

We may now invert this expression to find, 


— exp^‘V^(|x[ + Jjt'Dj + exp (i-Jis^x - x'|)) 


( 33 ) 


1 


G(x,x',t) = 7 = 

v ’ 2 


2 ia 

\x'\ + 2 ia 


( /, , , .,, 2 ^ 


exp i 


.(M+M) 


At 


+ exp 


f ( ^ 

(x-x) 


At 


JA 


(34) 


Solution for c with an Exponential Dependence on Time 

In this case we choose c(t)=-co for t<0 and c(t)=-co +(3exp(-at) for t>0, where c 0 is 
real such that c 0 >0, (3 is real, and ot>0. We choose the bound state for c(t)=-co as the 
initial condition where y/(x,0) = exp(-c 0 |jc[) . This represents a type of impulse. One 

can imagine the constant attractive bound state case, where at t=0, the potential jumps to - 
cO+b, then quickly falls back to -cO. This could represent a model for stimulated 
emission. 

We use the shift formula for Laplace transforms, which is given by, 


L[exp(at)f(t),s] = L[f(t),s - a ] 
By Equations 35 and 6 we find, 


(35) 


W{0,s) 


= J\§)__ _ftP y(0, J+g ) 

S~IC 0 'SIS-'VICq 


(36) 


To obtain a solution, we could iterate this expression repeatedly to obtain the series 
solution, 

(37) 


n-0 
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But since we already know the form of the series, we may substitute Equation 37 
directly into Equation 36. After doing this and solving for the a's by equating terms with 
like powers of (3, we find a n (s) is given by 


a n( s ) = 

a n{s) = 


o 


s-ic, 


o 


s + na - icg j=0 -Js+ joc - Vic 0 


,n >0 


(38) 


Discussion 

It was shown that the time dependent Schrodinger equation with a time dependent 
delta function potential can be solved exactly in many special cases. We were able to 
show that, in contrast to the case of sudden change, the adiabatic change in the potential 
leads to a conservation in the normalization of the probability density of the bound state. 
This is a rather significant result because it proves the system is reversible to within a 
phase factor if the change is slow enough. 

The periodic case and multi-delta function case was excluded for the sake of brevity 
and because they deserve a separate treatment due to the complexity of the solutions. The 
case of the sinusoidal time dependent potential can be treated in a similar manner as the 
exponential time dependent potential by representing the sine function as a sum of 
complex exponentials. Those two cases have also been treated elsewhere in a different 
context. As stated earlier the scattering problem was solved for the delta function with a 
strength that varied periodically in time [9] in as much as they were able to derive a 
transmission coefficient. We treated that same case as an initial value problem a decade 
before [8]. In the case of the multi-delta-function potential, this has also been recently 
treated in the literature [14,15] in the case of constant in time delta function potentials 
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and also in the case of delta function potentials moving with constant speed [13]. We also 
treated that case more than a decade before [8] for the case where the strength varies in 
time. 
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